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We present general arguments and construct a stress tensor operator for finite lattice spin models. 
The average value of this operator gives the Casimir force of the system close to the bulk critical 
temperature Tc- We verify our arguments via exact results for the force in the two-dimensional 
Ising model, d-dimensional Gaussian and mean spherical model with 2 < d < 4. On the basis of 
these exact results and by Monte Carlo simulations for three-dimensional Ising, XY and Heisenberg 
models we demonstrate that the standard deviation of the Casimir force Fc in a slab geometry 
confining a critical substance in-between is fci,r_D(T)(A/a'*~^)^/'^, where A is the surface area of the 
plates, a is the lattice spacing and D{T) is a slowly varying nonuniversal function of the temperature 
T. The numerical calculations demonstrate that at the critical temperature Tc the force possesses 
a Gaussian distribution centered at the mean value of the force < Fc >= kt,Tc{d — l)/S./{L/a)'^ , 
where L is the distance between the plates and A is the (universal) Casimir amplitude. 
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I. INTRODUCTION 

If material bodies are immersed in a fluctuating 
medium the surfaces of these bodies impose boundary 
conditions that select a certain mode spectrum for the 
fluctuations. This leads to a contribution into the ground 
state energy of a quantum mechanical system, or to the 
free energy of a critical statistical mechanical system, 
which depends on the geometrical parameters character- 
izing the mutual position of the bodies and their shape. 
\ This is known as the Casimir effect 0, 0, 13 • 

According to our present understanding, the Casimir 
effect is a phenomenon common to all systems character- 
ized by fluctuating quantities on which external bound- 
ary conditions are imposed. Casimir forces arise from an 
interaction between distant portions of the system medi- 
\ ated by fluctuations. 

In quantum mechanics one usually considers fluctua- 
tions of the electromagnetic field. In this case correla- 
tions of the fluctuations are mediated by photons - mass- 
less excitations of the electromagnetic field Q, |^ ^ |^ . 
In statistical mechanics the massless excitations can be 
generated by critical fluctuations of the order parameter 
around the critical temperature Tc of the system 0, ■ 
Goldstone modes (or spin wave excitations) in 0{n) mod- 
els at temperatures below Tc also provide massless exci- 
tations 0, 0, Fluctuations of this type are scale 
invariant and therefore the Casimir force is long ranged 
in the above cases. 

In this article we discuss the behavior of the thermo- 
dynamic Casimir force in systems with short-ranged in- 
teractions undergoing a second order phase transition. 

To be more specific, let us consider a statistical me- 
chanical system, a magnet or a fiuid, with the slab ge- 
ometry Tn~^ X L_\_, where d is the dimensionality of the 



system and periodic boundary conditions are applied. In 
the limit Ly — > 00 {L± fixed) the Casimir force per unit 
area is defined as 

/3Fcasi„.ir(T,T^) = _^M£iM^ (1.1) 

where fcx{T, Lj_) is the excess free energy 

/ex(T, Lj_) = f{T, Lj_) - Ti/buik(T) (1.2) 

of the system. Here f{T,L±) is the full free energy per 
unit area measured in units of ksT and /buik(T) is the 
corresponding bulk free energy density. 

According to the definition given by Eq. Hl.l|l the ther- 
modynamic Casimir force is a generalized force conjugate 
to the distance L± between the boundaries of the system 
with the property Tcasimir (T, L±) for L± — > 00. We 
are interested in the behavior of Tcasimir when Lj_ 3> a, 
where a is a typical microscopic length scale. In this limit 
finite size scaling theory is applicable. Then one has ^ 

(T, L^) = TI'*Xcasimir(T/Coo), (1-3) 

where ^oo is the true bulk correlation length, while 
Xcasimir IS an Universal scaling function. 

At the critical point Tc of the bulk system one has 
^oo = oo, and 

PcFc—iT,, L^) ^{d- 1) (1.4) 

where A is the so-called Casimir amplitude. This ampli- 
tude is universal, i.e. A depends only on the universality 
class of the corresponding bulk system and the type of 
boundary conditions used across L±_ . Obviously, one has 

A=Xcasimir(0)/(d-l). 
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The Casimir force may also be viewed from the point 
of view of conformal invariance of, e.g., critical systems 
|12| |. The Casimir force in its simple form for the film 
geometry (Ly oo, L± finite) is due to the L± depen- 
dence of the free energy f{T, L±) per unit area. The free 
energy therefore responds to any coordinate transforma- 
tion which changes the value of L±. On the other hand 
any coordinate transformation which transforms a slab 
of thickness L± into a slab of a different thickness is non- 
conformal. Therefore, the Casimir force is the response of 
the free energy f{T, L±) of the original slab to a noncon- 
formal coordinate transformation. The change of the free 
energy due to a nonconformal coordinate transformation 
is determined by the thermal average of the stress tensor 
tap associated with the Hamiltonian of the system [T5|. 
If the coordinate perpendicular to the surfaces of the slab 
is denoted by z, it is easy to prove, that the Casimir force 
in the slab is given by the thermal average of the stress 
tensor component tzz 113 (see below). 

Normally, one is interested in < t^z >, which de- 
termines the (average) value of the Casimir force. In 
addition, one can consider any realization tzz of the 
stress tensor to be proportional to any realization, e.g. 
instantaneous value of the Casimir force Fq-, where 
< Fc >= Fcasimir- That is the approach undertaken 
recently by Bartolo, Ajdari, Fournier and Golestenian 
[T^ . They consider a statistical- mechanical model of a 
d-dimensional medium described by a scalar field <& with 
an elastic energy density proportional to (V$)^, i.e. one 
considers an elastic Hamiltonian of the form 



V$(R) 



(1.5) 



where R = (r, z). They assume that the plates impose 
Dirichlet boundary conditions, i.e. $(r, 0) = ^(r, Lj_) = 
0. 

For the total Casimir force F^ = ijj^^Fc it has been 
found that: 
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F, 



Casimir 
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d-l 



cid) 



(1.6) 



where c(d) = (d - l)r(d/2)C((i)/(47r)'*/2 depends only on 
the spacial dimension d, while the variance of the force 
(AFq)'^, that can be considered of being produced of 
Na = {L\\/a)'^~^ independent strings is 



1 ^LnV-' 1 A 



/32 a-i-i ■ 



(1.7) 



In the above expressions A = L|| ^ is the cross-section of 
the system. From Eqs. (jl.fi|l and (|1.7|l one obtains that 
the " noise-over-signal ratio" p is 



<Fh> Uii 



(d-l)/2 



(d+l)/2 



(1.8) 



The probability distribution of the force has been found 
to be Gaussian, i.e. 



V{Fh 



x) 
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c) 



(1.9) 



The structure of the current article is as follows. First, 
in Section^we present some general arguments and con- 
struct a stress tensor operator for finite lattice spin mod- 
els. Then, in Section lTTll we verify our arguments by pre- 
senting exact results for the two-dimensional Ising model 
(see Section IIII All , d-dimensional Gaussian model (see 
Section IIII C|) , and for the mean spherical model with 

Monte Carlo results, based 



2 < d < 4 (see Section UlIBll 
on our definition, are given in Scction llVl for the behavior 
of the force and its variance. There the three-dimensional 
Ising f Section II VA|I . XY (Section IIV Bp and Heisenberg 
(Section llV ("P models have been considered. The article 
closes with a Discussion given in Section The set of 
technical details needed in the main text is organized in 
a series of Appendixes. 



II. THE STRESS TENSOR FOR LATTICE SPIN 
MODELS 

We will now reconsider the Casimir force for a d- 
dimcnsional anisotropic lattice 0{N) spin system with 
the Hamiltonian (see also Ref . [T^ ) 



(2.1) 



R fe=l 



where a d-dimensional simple hypercubic lattice with 
L||~^ X L_L lattice sites and Ly ^ Lj_ is assumed. The 
vector R indicates a lattice site and the vectors e^, 
k — 1,2, ... ,d connect nearest neighbor lattice sites on 
the simple hypercubic lattice. The spins Srt are con- 
sidered to be of 0{N) type. FoUowingRef . jl3l| and the 
general idea of conformal field theory |[lj|, we define the 
coupling constants in Ea. H2.1(l by 



Jfe(A) = J||(A) = J(e^), k^l. 



1, 



(2.2) 



which means that A = marks the isotropic point of 
Eq.lEU) due to Jy(0) = Jj_(0) = J(l). The function 
J{x) in Ea. H2.2|l is supposed to be smooth and monotonic 
in the vicinity of a; = 1 but is otherwise arbitrary. The 
critical point of the bulk spin model defined by Ea. (|2.1l) 
is given by an implicit equation of the type 



/C(/3cJi(A),...,/3,Jd_i(A),/3cJd(A)) = 
/C(/3eJ||(A),...,/3,J||(A),/3,Ji(A)) = 1, 



(2.3) 



where /?c = ^/iksTc). The function /C(ui, . . . ,Ud-i,u^) 
is a smooth function of d variables ui, . . . , Ud. Further- 
more, at /C(ui, . . . ,Ud) = 1 the function JC is invariant 
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with respect to any permutation of its arguments, be- 
cause the location of the critical point is independent 
of the labelling of the lattice axes. This also implies, 
that the derivatives dIC/duk, fc = 1, . . . , d for /C = 1 
all have the same value K,' ^ at the isotropic point 
Ml = U2 = ■ ■ ■ = Ud- For the two dimensional Ising 
model the function K is rigorously known 



/C(ui,U2) = sinh(2ui) sinh(2M2), 



(2.4) 



but for d > 3 exact results for K, are extremely rare. From 
Ea. (|2.3(l one immediately concludes, that in general the 
critical temperature Tc will depend on the anisotropy pa- 
rameter A. However, for the particular parameterization 
given by Ea. H2.2|l one finds at the critical point (see also 
Ref.0) 



= 



dK 



A=0 



= /3^(0)/C' [(d-l)J||(0) + Ji(0)] 

+ /3e(0)/C' [(d-l)J|^(0) + Ji(0)](2.5) 
= /3^(0)/C'dJ(l), 

which immediately yields /3c (0) = 0, i.e., for an infinites- 
imal anisotropy (A <C 1) the value of the critical temper- 
ature remains unchanged and is given by its value for the 
isotropic model. 

The correlation length of the anisotropic bulk spin 
system defined by Eas. H2.1| ) and (|2.2|l will also be 
anisotropic. In the vicinity of the (isotropic) bulk critical 
temperature one finds 

e||(A,0 =e||.o(A)Kr'' and e±(A, i) - a,o(A)|^r^ 

(2.6) 

where t = iT — Tc)/Tc is the reduced temperature and 
V is the correlation length exponent which is universal, 
i.e., independent of the anisotropy parameter A. At the 
isotropic point A = one has ^i[,o(0) = '?_L,o(0) = ^o- To 
simplify the notation we ignore the fact that the correla- 
tion length amplitudes in general also depend on the sign 
of the reduced temperature t. We therefore assume that 
t > in the following. 

In order to be able to apply finite-size scaling in the 
critical regime with respect to a single correlation length, 
say, , we employ the following anisotropic rescaling of 
the spatial coordinates: 



x'y.=Xk, k = 1, . . . ,d — 1, and 



eii,o(A) 
e±,o(A) 



(2.7) 

This transformation has the desired property, namely 
^11 ~ and = as can be easily verified from 
Eas. H2.6|) and H2.7|l . For the lattice sizes Ly and L± we 
find accordingly 



Lu = L\\ and L', 



eil.o(A) 
e±,o(A) 



Li 



i?(A) L_L. (2.8) 



The link between the explicit A-dependence of the free 
energy of our spin system according to Eas. (|2.1| ) and (|2.2|) 



and the Casimir force defined by Ea. Hl.l() is provided by 
Ea. (|2.8|) . In the limit Ly ^ oo we find (see also Eqs.(24) 
and (33) of Ref. O) 



Cjfcx 



lim 



I3J'{1) 



A=0 



(2.9) 



.fe=l 



where (...) denotes the thermal average with respect to 
the Hamiltonian given by Ea. H2.1|l at the isotropic point 
A = 0. From Eas. ljl.ip and 1)2. 9|l and the finite size scaling 
form 



9c 



t [LjU^of"] (2.10) 



of the excess free energy fex(t, L±) in the limit L|| — > oo, 
we obtain an expression of the Casimir force which is 
derived in detail in Appendices IXI and IbI From Ea. (|B8|) 
derived in Appendix ^ the operator form of the stress 
tensor component tj_j_ (R) can be read off as 



iii(R) = f3J\l)[R'iO)]-' 

-d-l 

X ^ SnSn+ek ~ {d — l)S'RS'R+ed 



.k=l 



(2.11) 



where Eq. (|B9|) was used and the operators Ti and Tit, are 
properly normalized Hamiltonians 7i(0) fsee Ea. (|2.1|l and 
Appendix 0) . 

Ea. H2.11() provides the connection between the stress 
tensor component t±± parallel to the surfaces of the 
slab and the spin lattice model given by Ea. H2.1|l . It 
is valid also for temperatures T > Tc and thus general- 
izes Eq.(36) of Ref.p'l. For T < Eq.lEHJl holds also 
for 0{N) symmetric spin models, because the correla- 
tion length ratio remains finite at T = and it 
can be continued analytically into the Goldstonc regime, 
where it can be used for the anisotropic rescaling of the 
coordinates according to Ea. H2.7|l . Finally, we note that 
Eg. (|2. 11(1 only holds for periodic boundary conditions. 

For the purposes of this investigation Ea. H2.11() serves 
as a prescription to obtain the universal scaling func- 
tion of the Casimir force in critical slabs with periodic 
boundary conditions by Monte-Carlo simulations. How- 
ever, Ea. H2.11l) contains the ratio ^||,o(A)/^_l,o(A) of the 
correlation length amplitudes for Ea. H2.1|l as a prefactor. 
For the two dimensional Ising model this function is given 
by [13 



eii.o(A) _ J^(A) + J|[(A)sinh[2/3e(A)J^(A)] 
U.oW J\\ (A) + J±(A) sinh[2/3,(A) J|| (A)] 



(2.12) 



and by virtue of Ea. (|2.2|l Eg. (|2. 11(1 can be made explicit. 
But in d > 3 no such information is available. However, 
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in the critical regime, where Ea. (|2.11|l will be applied, the 
scaling argument will typically not exceed values 

of the order 10. This means that we will be dealing with 
reduced temperatures in the range |t| < 10(ij^/^o)~^^''i 
i.e., the relevant temperature range diminishes as the 
system size Lj_ increases. As can be seen explicitly in 
Eq. (|2.12|l the correlation length amplitude ratio only de- 
pends on the temperature T and does not display any 
scaling behavior. Therefore, the generally unknown pref- 
actor in Ep. H2.lll) can be treated as a constant for suffi- 
ciently large system sizes which can be determined by a 
normalization of {t_\__\_) to known results aXT — Tc 6]. 

We end this section with some observations that turn 
out to be very helpful for analytical calculations of the 
variance of the force. Since they are model indepen- 
dent we give them before passing to explicit calcula- 
tions presented in the next section. Lets us consider an 
anisotropic Hamiltonian of the type given by Eq. H2.1|l . 
where Ji(A) = • • • = Jd_i(A) = J|| (A) = (1 + A) J and 
Jd{X) ~ J^-{^) = ^ {d— 1)A) J. Then, it is easy to see 
that 



H(A) =H(0)-AJ^ 



R 



or, equivalently, 
where 



(2.13) 



,fc=i 



A^f_^,^(R), 

R 



(2.14) 



R+et - (d - l)S'R5R+ed 



.fc=l 



(2.15) 

differs only by a multiplying factor from the stress tensor 
ij^,_L(R) defined in Eq. ifTTH) . 

Let /(T, A) is the total free energy per unit spin of a 
system with the Hamiltonian (|2.14|l . Then, taking into 
account the translational invariance symmetry, it is easy 
to see that 



«x,x(R)> ^ 4 

= l3J{d 



Pf{T, A) 
- 1) [{SoS,, 



(2.16) 



A=0 



Therefore, in order to calculate the average value of 
'^-L,-l(R-) one needs either to know the finite-size free 
energy density of an anisotropic system, or, what is 
much simpler, the nearest neighbor two-point correla- 
tions along the axes of the isotropic finite system. 



III. ANALYTICAL RESULTS 

In this section we summarize our analytical results for 
the two-dimensional Ising, for the spherical model with 
2 < d < 4, and for the Gaussian model. Their derivation 
for the Ising model is given in Appendix El while ones 
for the spherical model are given in Appendix IdI 



A. Two-Dimesnional Ising Model 

For the two-dimensional Ising model on a square lattice 
with geometry L x M the lattice representation of the 
stress tensor is well known for a long time 
our notations, using Eqs. (|2.11() and H2.12|l . we obtain 



/3J'(1) 



\- 






dX 




X=0- 



f3J{l + sinh(2/3J)) 



2(sinh(2/3J) - 2/3Jcosh(2/3J) - 1) 



(3.1) 



At the critical point (3c of the isotropic system one has 

El 



1 = sinh(2/3,J), 



(3.2) 



and the right-hand side of the above equation simpli- 
fies essentially becoming simply — 1/(2\/2). Therefore, 
at T = Tr, the stress tensor reads 



^x^xi^j j) 



1 



2V2 



(3.3) 



which is exactly the form considered in . In the limit 
M ^ 00 at the critical point Tc of the bulk system one 
has fil 



6 



(3.4) 



where c = 1/2 is the so-called central charge of the Ising 
model The Casimir amplitude is 12J 



1 

A = -- c, c = -. 
6 ' 2 



(3.5) 



It is easy to see that close to Tc the right-hand side of 
Eq. (|3.1|) becomes 

il + ^^)+0{iP-Pcr). (3.6) 



2V2 



Since i/ = 1 for 2d Ising model, it is clear from Eq. H3.4(l 
that the contributions to the Casimir force due to the 
term proportional to P—Pc m the above expression will be 
of the order of L~^. Such contributions will be neglected. 
Therefore, in the critical region of the finite system we 
conclude, that the stress tensor is given by 



tx,x (^7 j ) 



1 



2V2 
1 



iPc-P)in-Hb). 



(3.7) 



One can interpret the variance of the stress tensor 
'^tx.xihj) as a variance of a local measurement of the 
Casimir force made near the point (i,j). For the lead- 
ing behavior of the variance at Tc one then has (see Eq. 

m) 



(3.8) 
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Definitely, in addition from the above nonuniversal part 
the variance contains also universal parts that are negli- 
gible in comparison with the nonuniversal one. 

As we said above, we will interpret {tx,x{i,j)) as a 
local measurement of the Casimir force made near the 
point Let us imagine, that we are collecting mea- 

surements from all the points belonging to the " surface" 
(l; j)i J G [1; ■ ■ ■ I (in the very same way one can con- 
sider the opposite "surface" {L,j), j e [1, • ■ ■ ,M]). The 
surfaces are important because they are the only place 
where in an experiment the Casimir force is experimen- 
tally accessible. To characterize the force measured on 
the whole surface, instead oftx^xih j), one has to consider 
X]j tx,x{^,j)- Taking, in a first approximation, any local 
measurement to be independent from the other ones, one 
obtains that 

Aj^txA^j) - MAt,,,(l,l) ~ 0.363 M, (3.9) 

j 

which implies that, indeed, in agreement with p^ . 

{A/3F^f (X N^-^ ^ {Li\/af-\wheTe d ^ 2. (3.10) 

An estimation can be also derived for A'^^ ^tx^x{i, j)- 
With a variance of such a type one deals when, say, Monte 
Carlo simulations of the force are performed. One obtains 
(see Eq. HCWi ) 



A^tx,x{'i-'j) 



ML ~ 0.068ML. (3.11) 



We again observe that the variance of the sum of tx^xihj) 
is proportional to the total number of summands in 
this sum. The coefficient of proportionality for 2d Ising 
model, when the sum is over all spins in the finite system, 
turns out to be 0.068. 

Unfortunately, as far as we are aware, the finite-size 
properties of the free energy of the two-dimensional 
anisotropic Ising model under periodic boundary condi- 
tions are not available for T ^Tc- This makes the com- 
parison of the direct derivation of the force as a derivative 
of the finite-size scaling excess free energy and as average 
of the operator H2.11|l a challenging task. Even more - 
the behavior of the finite-size free energy of the isotropic 
system is o nly known for moderate values of the scaling 
arguments T?!. Nevertheless, from [l^ one can extract 
the following results for the scaling functions of the excess 
free energy and the Casimir force 

• excess free energy 
The scaling function of the excess free energy is 



i=2 ^ ' 

where — tt < x < 2tt, and the scaling variable is 

X = SKrtL. 



• Casimir force 

The scaling function of the Casimir force Xcasimir 
is related to that one of the excess free energy via 

d 

^Casimir = -'^ox(a;) - X — Xcxix). (3.13) 

OX 

Then, from H3.12|l . one immediately obtains 



a 



1 9 ^ 



1=2 



1/2 



(3.14) 



(|-)''(l-2z)(l-2-2'+i)C(2z-l). 



B. The Spherical Model 

We consider a spherical model on a c?-dimensional hy- 
percubic lattice A G Z"^, where A = Li x L2 x ■ ■ ■ Ld- 
Let Li = Nitti, i — 1, • ■ • , d, where Ni is the number of 
spins and is the lattice constant along the axis with 

being a unit vector along that axis. With each lattice 
site r one associates a real- valued spin variable Sr which 
obey the constraint 



(^2) = N, 



(3.15) 



where N = N1N2 ■ ■ ■ Nd is the total number of spins in 
the system. The average in H3.15II is with respect to the 
Hamiltonian of the model which is 



m = -|/3y 5r J(r, r')Sr' + S ^ 



(3.16) 



In the current article we will consider only the case of 
nearest-neighbor interactions, i.e. we take J(r, r') = 
J(|r — r')| = Jj, if r — r' = ±ej, i ~ 1, • • • ,d, and 
J(r,r') = otherwise. 

For such a model it can be shown that, under pe- 
riodic boundary conditions, the free energy of the model 
(per unit spin) is given by 



f3fiK,N) = 



K 

In K 

2tx 



+ — y In 
2iV ^ 



sup 

io>0 



1 



1 - 



2 

m 

J(0) 



Kw (3.17) 



while the two-point correlation function is 



ui + 1 - J(k)/J(0) 



(3.18) 



Here s = K{w + l)/2, K = /3J(0), where J(k) is given 
by the Fourier transform of the interaction 



J(k)=yj(r)e 



ik.r 



(3.19) 
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and the wave vector k = {fci, A:2, • ■ ■ , fcrf} G ^Ba is with 
components fc^ — 27mi/Li, where rii = , Ni, i = 

1 , • • • ,d. The equation for the spherical field w reads 



N ^ 



1 



w + 1 - J(k)/ J(0) 



(3.20) 



which leads immediately to G(0,f, N) = 1. 

For nearest neighbor anisotropic interactions it can be 
shown that (see Eq. (jDlGp ) 



(3.21) 



where f j is the correlation length in direction j, and hj = 
"^jj which leads to the following explicit form 
of the stress tensor within the spherical model (see Eq. 

(Imsli 'i 



1 



.k=l 



+-^{.Pc-P)[n-Hb 

a V ^ 



(3.22) 



Here Ti. is the Hamiltonian (normalized per unit particle) 
of the finite system, and Ti.b is that one of the infinite 
system. 



1. Evaluation of the finite-size excess free energy of the 
anisotropic system 

First, one can demonstrate that the critical coupling 
of the anisotropic bulk system is (see Eq. l)D30|l 1 

d poo ^ 

K, = 2/3, ^ J, = W^d(0|b) = / dx\{ e-"^^/o(a;6j). 

(3.23) 

Then, close to K — Kc, when 2 < c? < 4, for the scaling 
form of the excess free energy /3(/ — fb) (per spin) in 
the limit of a film geometry Ni,N2, - ■ ■ , Nd-i oo one 
obtains (see Eq. ID35I) 

P[f{K,N^\h)~fb{K\h)] = |ixi(y-yoo) (3.24) 



i r(-d/2) /fo^Y''"'^/" / 

2 f47r)'i/2 Uii j V 



{d-l)/2 oo 



2 



(27r)''/2 V^>j| 
In the above equation 



„d/2 



9=1 



rf- 2' 



(3.25) 



is the temperature scaling variable, j/oo = '2wbN'j_/b±^ is 
the solution of the bulk spherical field equation (see Eq. 

(jmU) 

-^-x,-^S^^^('A'"'''\'''-' (3 26) 
while ?/ = 2wN'^/b± is the solution of the finite-size 



spherical field equation (see Eq. (jDSSj) 

1 r(i-d/2) fb^Y"-'^/' ,^,^, 



2^' 



(27r)'i/2 \ b 



b 

(d-l)/2 



y 



r'- - (3.27) 

d/4-l/2Y- ^<i/2-l(9%/y) 



9=1 



d/2-l 



For the Casimir force one derives (see Eq. (|D40|l ) 
1 



asimir -^_L 



-d 



7Xl{y-yoo) 



(3.28) 



(d-l)/2 



i r{-d/2) 

2 (47r)''/2 



d/2 



d/4 



E 

9=1 



Kd/2{qy/y) 



^d/2 



whereas for the Casmir amplitudes A we derive (see Eq. 
lin44ll ^ 



d(27r)rf/2 



V2^ 



Kd/2+i{q^/W) 



9=1 



,(i/2-l 



(3.29) 



with /3c-Fcasimir(J^c,i_L) = (rf - 1) A L^''. The exact 
value of Uc and A in an explicit form is only known for 
d = 3. Then 



2/c-41n2[(l 



(3.30) 



(this value is well-known and seems that has been derived 
for the first time in 20]), and, then, one obtains 



2C(3) 



(3.31) 



This is the only exactly know Casimir amplitude for a 
three dimensional system. 



2. Evaluation of the average value of the stress tensor 

For the scaling form of the average value of the stress 
tensor one derives (see Eq. (ID54II ') 



2{d-l) 
(27r)'i/2 



9=1 



jd/2 



L,(d+2)/4 Kd/2-l{qVy) 



E 

9=1 



„d/2-l 



+ ^^^^i{y-yoo)\Nz''. (3.32) 
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It is also possible to demonstrate (see Appendix ID|| that 
the above expression is equivalent to /3Fcasimir given by 
Eq. (|3.28|l for the isotropic system (when fey = b±), i.e. 
indeed, 



(tj_j_(R)) — /3i^Casimir 

for the spherical model. 



(3.33) 



Evaluation of the variance of the stress tensor 



we conclude that in the Gaussian model the stress tensor 
again is 



(3J 



d-l 



.fc=l 



-^iPc-P){n-nb), (3.37) 



where Ti. is the Hamiltonian (normalized per unit parti- 
cle) of the finite and Hb of the infinite model. 



For the variance of the Casimir force in the spherical 
model at T = Tc one obtains (see Eq. H3.34|l ') 



At I 



A J2 i^.^(R) 
ReA 



0.107 iV_LA^|f 



(3.34) 



This will be also the leading result everywhere in the crit- 
ical region - since it is coming from the nonsingular part 
of the free energy around T = an analytical expansion 
should be possible. It is also clear that if the summation 
over R in 13.34|1 is not over the total number of particles 
in A, which is N^N^"'^, but over, say, all the spins from 
one of the boundary, then the corresponding variance will 
be proportional to the total number of spins there, i.e. 
that 



(A/3F^)2 cx TV, 



(3.35) 



exactly as it has been found in 14], see Eq. 11.7|l . 



C. The Gaussian Model 

In order to simplify the notations we define the Gaus- 
sian model in the same way that we have defined the 
spherical model. Actually the spherical model is a Gaus- 
sian type model with one additional constraint, given by 
Eq. H3.15() fixing the average length of all the spins in 
the system. To be more precise, we suppose that the 
Hamiltonian H of the model is again given by Eq. I|3.16|l . 
where, as before, s = K{1 + w) /2, an d K = /3J(0). The 
only difference is that now Eq. H3.15|l is missing and w is 
not a quantity which behavior has to be derived from it, 
but a parameter which describes the deviation from the 
critical point, i.e. w = (/3c — P)/f3- As a result, the free 
energy density of the model becomes 



/3/(if,N) = i 



27r 



C/(w,N) - ^Kw, (3.36) 



where U{w, N) is given by Eq. I|D3|I . while the two-point 
correlation function is still determined by Eq. H3.18|l . 
Then, for a system with anisotropic short-range interac- 
tion, proceeding in the same way as in Section IlII Bl we 
derive that Eqs. (|D13|I - ljD17|l are still valid, wherefrom 



1. Evaluation of the finite-size excess free energy of the 
anisotropic system 

The analysis of the excess free energy cam be per- 
formed along the same lines as in the case of a spherical 
mod el. Fo r example, for U{w,N) the Eqs. (|DT9)| - ||D23)) 
and (|D28|I are still valid. On the basis of these equations 
we immediately obtain that in the case of a Gaussian 
model the excess free energy in a film geometry is 



p[f{K,N^\h)-MK\h)] 



(3.38) 



{d-l)/2 oc 



y 



9=1 



d/2 



In the above equation the temperature scaling variable is 



y = 2wN\_^'' /bj_, where v = 1/2, and w = Pc/P - 1- 

(3.39) 

From Eg. H3.38|l . using the property of the K,y{x) func- 
tions |l3 that 

^ [y''K,,{ay)] - -ay'' K,^,iay), (3.40) 

we immediately derive that the Casimir force in such an 
anisotropic Gaussian model is given by 



x(d-l)/2 

(y) (3.41) 



(27r)'i/2 



9=1 



_^^rf/4+l/2^ 



qd/2-i j 



Note that despite the similarities with the spherical 
model both the excess free energy and the Casimir force 
differs essentially for the two models. Let us demonstrate 
that even more explicitly on the example of the Casimir 
amplitudes A. We remind that in he spherical model 
they are given, for 2 < d < 4, by Eq. (|3.29|l where yc 
is the solution of the spherical field equation at P = Pc- 
In explicit form we have been able to solve this equation 
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and to calculate A only for d = 3. The situation with 
the Gaussian model is much simpler. At the critical point 
f3 = Pc one has y = 0, and, therefore, from Eq. H3.41(l . 
or Eq. (|3.38|) . we obtain (in the isotropic system) 



A 



r{d/2)c{d) 

_d/2 



(3.42) 



So, for d = 3 one has A = —({3)/2tt and, therefore. 



^Spherical Model AQauggian Model: 



d = 3. 



(3.43) 



2. Evaluation of the average value of the stress tensor 

Having in mind that u = -§p{fif) and using 1)3.40(1 . for 
the difference of the finite-size and bulk internal energy 
densities one can easily derive from (|3.3t)|) and 



U — Ub 



MP y'^/^+^/^g^ j^./2-i(gV^) ,,_, 



(27r)'i/2 (/3, - /3) ^ g-^/^-i 



(3 44 ) 

where y = 2dwNl and w — Pc/P — 1- Next, from H3.18|l . 
or H3.36|l . for the stress tensor (|3.37|) of the Gaussian 
model we derive that 

(t (R)) - cos(fciai) - cos(fcdQrf) 



Pc/P 



./4+i/2g ^^/2-iy) ^-.^ (3.45) 



q=l 



where we have taken into account that v =^ 1/2. Apply- 
ing to the first row in the above equation the same way of 
acting, as in the case of the spherical model, and replac- 
ing Pc/P by 1 (since we are close to the critical point), 
we derive 



^ i yd./4+l/2 

d(2^)^/2 Y ^ 

q=l ^ ) 



Kd/2~iiq^/y) 



,d/2-l 



iV7 



(3.46) 



Now it only remains to show that the right-hand side of 
the above equation is indeed equal to the right-hand side 
of Eq. H3.41II (for b± — b\\), which gives the Casimir force 
calculated in a direct manner as a derivative of the finite- 
size free energy with respect of the size of the system. In 
order to demonstrate this, let us note that, according to 
the identity (|D43|) . 



(3.47) 



Kd/2+l{x) Kd/2-l{x) + -Kd/2ix). 



Inserting 13.47|l in H3.46|l and comparing the result with 
Eq. ((3.41(1 . we conclude that, indeed. 



(tj_j_(R)) — /3-Fcasimir 

for the Gaussian model. 



(3.48) 



3. Evaluation of the variance of the stress tensor 

For the variance of the stress tensor all the equations 
from the corresponding part of Appendix^for the spher- 
ical model are still valid. That is because the leading 
contribution of the variance is stemming from the regu- 
lar part of the bulk free energy Ub{0\h) (see Eq. I(D25() ') 
evaluated at T = Tc (see Eq. ((D65|l ). This observation 
leads to the conclusion that, as in the spherical model 



A<j_,j_ = A ^ i±,±(R) ^ 0.107 iVj.iV|, (3.49) 
rga 

where the summation is over all the spins of the system. 



IV. MONTE-CARLO RESULTS 

The foundation of our Monte-Carlo investigations of 
the critical Casimir force is laid by Eos. ((2. 11(1 and ((lj8(l . 
respectively. Apart from the a priory unknown coefficient 
J'(l)/i?'(0) and the bulk energy density Ub the numerical 
evaluation of Eq. ((2.11(1 is absolutely straightforward and 
apart from usual algorithmic precautions in the critical 
regime no special techniques are required. However, as 
has become obvious in, e.g.. Eg. ((3. 11(1 . the statistical er- 
ror of the estimate will increase with the system size if 
the number of Monte-Carlo sweeps is kept constant. In 
order to approach the asymptotic regime larger system 
sizes, say, TVy = 120 and N± = 20 lattice sites in d = 3 
are required which means that a reliable estimation of the 
Casimir force remains computationally demanding as far 
as the required CPU time is concerned. 

We employ a hybrid algorithm '22] which consists of 
Metropolis (231] and single cluster updates for the 
Ising model, for XY and Heisenberg simulations over- 
relaxation updates j25j are employed as a third update 
method. Cluster updates are only used in the immediate 
vicinity of the critical point, e.g., for —0.02 < t < 0.02 
for the system size indicated above. Typically, we have 
performed between 4.8 x 10^ and 9.6 x 10^ Monte-Carlo 
steps per spin. In order to cope with the high demand 
of CPU time for larger systems we have performed part 
of our simulation in parallel on two-processor Intel Xeon 
system and on a four-processor DEC Alpha system using 
the OpenMP Standard for SMP programming. A few 
runs have also been performed on a two-processor AMD 
Optcron system. 

We first investigated the energy dependent contribu- 
tion (/3c (0) - P){u~ Ub) to Eqs.(Eni and (jBSl by a se- 
ries of simulations on a cubic geometry for A^n N± = 
20 ... 80 in order to obtain reliable estimates for the 
bulk energy density Ub- It turns out, that within the 
range —0.2 < t < 0.2 of reduced temperatures, vari- 
ous aspect ratios N\\/N± = 3,4,6,8, and several system 
sizes A'^ii = 60 . . . 120 the energy dependent contribution 
(/3c(0) — P){u — Ub)/dv is always negligible. As a typical 
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result we obtained that for forces of the order 10""'^ with 
a statistical error in the range 10^^ . . . 10"'^ the energy 
contribution (/3c(0) — f3){u— U}j)/dv remains in the range 
10-4 . . . 10-5 for all models. The prefactor J'(l)/i?'(0) 
roughly evaluates to J'(l)/i?'(0) ~ 0.3 in all cases. We 
therefore conlude that we can safely ignore the energy 
contribution to the Casimir force for our Monte-Carlo 
investigations of the Ising, the XY, and the Heisenberg 
model in three dimensions. 



A. Three-Dimensional Ising Model 

As expounded above, we have neglected the en- 
ergy dependent contribution to Ea. (|B8|) for our Monte- 
Carlo evaluations of the scaling Function 9per{x), x ~ 
tiL^/a)^/" of the Cas imir force. From extended simula- 
tions for various aspect ratios N^\/N± = 3, 4, 6, and 8 we 
have arrived at the conclusion that corrections to 9per (x) 
due to finite aspect ratio are by far negligible within the 
statistical error for A^||/A^j^ = 6. In fact, our results for 
A'^ll /Nj_ — 4 can hardly be distinguished from correspond- 
ing results for larger aspect ratios. We have therefore 
fixed the aspect ratio to the value 6 and performed simu- 
lations for iVj^ — 16, 20, 24, and 30. The resulting scaling 
plot for Oper{x) is shown in Fig^ For T > Tc finite-size 




FIG. 1: Scaling function 9per{x) of the Casimir force for the 
d = 3 Ising model in a slab geometry for periodic bound- 
ary conditions as function of the scaling variable x = tN^^'^ , 
where A^x — L±/a is the number of lattice layers. The as- 
pect ratio is chosen as N^^/Nj_ = 6 (see main text). Finite 
size scaling according to our expectation is confirmed within 
two standard deviations, where u = 0.63 has been chosen. 
The vertical scale has been adjusted according to the esti- 
mate 6'per(0) = 2Aper,n=i = -0.306 (see Ref..26]). The error 
bars displayed here correspond to one standard deviation. 



scaling works very well, whereas for T < Tc data col- 
lapse for N± = 30 is not as good. However, the data 
collapse improves upon increasing the statistics and so 
we have performed 9.6 x 10^ Monte-Carlo steps per lat- 
tice site for the largest lattice iVj_ = 30 for T < T^. With 



the estimate v = 0.63 for the correlation length expo- 
nent we finally obtain scaling within two standard devia- 
tions. The scaling function 9per{x) decays exponentially 
for X —^ ±oo and has its minimum below Tc- However, 
due to the magnitude of the statistical error its loca- 
tion cannot be determined accurately enough to exclude 
X = with reasonable certainty. The Monte-Carlo data 
for the Casimir force are not normalized due to the a pri- 
ory unknown prefactor J'(l)/i?'(0) in Ea. HB8|l . The data 
displayed in Fig^ have therefore been scaled in such a 
way that 0per(O) = 2Aper,n=i is given by the best known 
-0.153 for the Casimir anmplitude 



peVjU—l 



mpl i 
(26|. 



estimate A 

Aper,n=i for the three-dimensional Ising model l2 

The scaling function displayed in FigQ] has been ob- 
tained from Ea. ljB8|l . where a spatial average over all lat- 
tice sites is performed. As expounded in Sec. Ill (see also 
Ref.jQ) this leads to a certain size dependence of the 
variance of the stress tensor as, e.g., given by Ea. H3.34|) 
for the spherical model and by Ea. H3.49|l for the Gaus- 
sian model. In order to investigate the variance also for 
the Ising model in d = 3 we have recorded the distri- 
bution function of the stress tensor during our Monte- 
Carlo simulations. It turns out that the shape of the 
distribution function is captured by a Gaussian distribu- 
tiuon to a very hig h degree of accuracy also for the Ising 
model (see Ref.|Ij|). We are therefore able to extract the 
variance of the stress tensor average from a least square 
fit of the measured distribution function to a Gaussian, 
where the variance is one of the fit parameters. Guided 
by Eas. H3.34|) and (13.491) we have normalized the vari- 



ance to Nj_ in order to obtain a linear law at fixed as- 
pect ratio. Our results for iV||/7V_L = 6 are displayed 
in Fig[21 The functional dependence is indeed linear 




FIG. 2: Variance At^.x of the stress tensor for the Ising 
model in d = 3 (see Eas. H3.34^ and H3.49|l '). normalized to 
A^x at fixed aspect ratio N\\/N± = 6 as function of A'"x for 
different reduced temperatures t in the critical regime. The 
behavior is linear as indicated by the straight lines connecting 
the data points. Their slopes have been evaluated as 1.39 for 
t = -0.01, 1.55 for t = 0, and 1.61 for t = 0.01. The statistical 
error of the data (one standard deviation) is smaller than the 
symbol size. 
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and the slope at t = {T = Tc) is 1.55 as compared 
to 0. 107(i V||/jVj _)^ ~ 3.85 for N^/N^^ = 6 according to 
Eg s. 1)3.34(1 and ((3.49(1 for the spherical and the Gaussian 
model. The strict linearity also prevails for other tem- 
peratures in the scaling regime as shown in Fig|2| The 
quadratic dependence of At±,±/N]_ on the aspect ratio 
has also been confirmed for the Ising model in c? = 3 at 
T = Tc from simulations at different aspect ratios (not 
shown) . 



B. Three-Dimensional XY Model 

In accordance with our findings for the Ising model we 
find that the value 6 for the aspect ratio of the simula- 
tion lattice is also a good choice for the XY model. We 
have performed simulations for iVj^ = 16, 20, 24, and 30, 
where 4.8 x 10^ Monte-Carlo steps per lattice site have 
been performed for all lattice sizes. It turns out that the 
energy dependent contribution to Ea. ((B8(l can again be 
disregarded within the statistical error obtained from the 
simulations. 

As in the Ising case we determine the normaliza- 
tion factor J'(l)/i?'(0) in Ea. ((B8l) from the requirement 
^per(O) — 2Aper,n=2- Unfortunately, all estimates for 
Aper,n=2, which are currently available, are based on 
the e-expansions quoted, e.g., in Ref.|^. Independent 
Monte-Carlo estimates for the Casimir amplitudes of the 
XY model do not exist and rigorous results for the two 
- dimensional XY model are limited to temperatures be- 
low the Kosterlitz - Thouless Temperature, where the 
model renormalizes towards the two - dimensional Gaus- 
sian fixed point. The Gaussian model in d = 2 is char- 
acterized by the central charge c = 1 and therefore the 
Casimir amplitude for the two - dimensional XY model 
in the low temperature limit is given by Aper,n=2.d=2 — 
-7rc/6 = -7r/6 ~ -0.5236 [23. 

Apparently, the e-expansion underestimates the mag- 
nitude of the Casimir amplitude Apcr,n=\ of the critical 
Ising model in d = 3, i.e., e = 1. From the structure of 
the critical Ginzburg - Landau 0^ theory and the nature 
of the two - loop approximation to the Casimir amplitude 
we expect that the e - expansion will also underestimate 
the magnitude of Aper,n for any n. This leads us to the 
conclusion that the e-expansion of the ratio 



(4,1) 



is more accurate in d = 3 than the £-expansion for nu- 
merator and denominator individually (see Ref.[^). We 
therefore adopt the approximation 

5 /n-F2 V ' 
+ 8 " 3. 



A, 



^0.153 n 



1 - 



(4.2) 



as our estimate for Aper,n in d 



A 



per,n— 1 



Ea. ((4.2() we then have 
A. 



3 in the following, where 
0.153 (see above) has been used. From 



for the three dimensional XY model. The resulting scal- 
ing plot for 9per{x) is shown in Fig|3| For T > Tc finite- 




FIG. 3: Scaling function dpsr{x) of the Casimir force for the 
d = 3 XY model in a slab geometry for periodic boundary 
conditions as function of the scaling variable x = tN^'^ , 
where A''_i = L±/a is the number of lattice layers. The as- 
pect ratio is chosen as N^^/N± = 6 (see main text). Finite 
size scaling according to our expectation is confirmed within 
two standard deviations, where v — 0.67 has been chosen. 
The vertical scale has been adjusted according to the esti- 



mate Sper(O) 



2A 



per, 71—2 



-0.56 (see main text). The error 



■per, 2 



-0.28 



(4.3) 



bars displayed here correspond to one standard deviation. 



size scaling works very well, whereas for T < Tc data 
collapse for A^j^ = 30 is again not as good. However, 
the data collapse ist still acceptable within two stan- 
dard deviations, so we refrain from performind additional 
runs here. The scaling function Operix) decays exponen- 
tially above for a; ^ 00 and displays a minimum be- 
low Tc- Unlike the Ising model the XY model exhibits 
lang-ranged correlations also below Tc (goldstone modes) 
which are a prominent feature in Fig|3| The scaling func- 
tion Operix — !■ —00) saturates at about half its minimum 
value and does no longer decay to zero. 



We have also evaluated the size dependence of the vari- 
ance of the stress tensor for the XY model along the lines 
of the previous analysis for the Ising model. The distri- 
bution function of the stress tensor is again given by a 
Gaussian to a very high accuracy. The corresponding re- 
sult for Atj_^j_/A^^ is shown in Fig^ The functional de- 
pendence is again linear and the slope at t — (T — Tc) 
is 2.78 as compared to 0. 107(A ^||/A_ l)^ - 3.85 (see previ- 
ous subsection and Eas. ((3.34(1 and ((3. 49(1 1 for the spher- 
ical and the Gaussian model. The strict linearity also 
prevails for other temperatures in the scaling regime as 
shown in Fig0] In summary the XY model behaves just 
as the Ising model with respect to the variance of the 
stress tensor. 
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FIG. 4: Variance At a 



in d = 3 (see Eas. l3.34|l and 113. 49H . normalized to A^^ 



of the stress tensor for the XY model 

a,t 

fixed aspect ratio Ni\/N± = 6 as function of N± for different 
reduced temperatures t in the critical regime. The behav- 
ior is linear as indicated by the straight lines connecting the 
data points. Their slopes have been evaluated as 2.63 for 
t = -0.01, 2.78 for t = 0, and 2.86 for t = 0.01. The statisti- 
cal error of the data (one standard deviation) is smaller than 
the symbol size. 



C. Three-Dimensional Heisenberg Model 

We have repeated the simulations finaUy for the 
Heisenberg model in c? = 3 with the same geometric and 
statistical data as for the XY model for the same reasons 
discussed above. We note again that the energy depen- 
dent contribution to Ea. (|B8(l can be disregarded within 
the statistical error obtained from the simulations. 

The normalization factor J'(l)/i?'(0) in Eq.(|B8l) is 

i\ A 

per,n— 3 ? 



determined from the requirement 
where the estimate 



.(0) 



2A, 



-0.39 



(4.4) 



used here has been obtained from Ea. l4.2|l for n = 3. 
The resulting scaling plot for 6per{x) is shown in Fig|Sl 
For T > Tc finite-size scaling works very well, whereas 
for T < Tc the scatter of the date is larger than for 
the XY model. However, the data collapse ist still ac- 
ceptable within two standard deviations. The qualitative 
shape of the scaling function 9per{x) is the same as for 
the XY model. The Heisenberg model also exhibits lang- 
ranged correlations below Tc (goldstone modes) which 
are a prominent feature in FigEl The scaling function 
Opcr{x — > — oo) saturates at about three quarters of its 
minimum value. 

Finally, have evaluated the size dependence of the vari- 
ance of the stress tensor for the Heisenberg model along 
the lines of the previous analyses for the Ising and the XY 
model. As before the distribution function of the stress 
tensor is given by a Gaussian to a very high accuracy . 
The corresponding result for At±^±/N^ is shown in FigEl 
The functional dependence is linear and the slope at t = 
(T = Tc) is 3.92 as compared to 0.107{NJNj_)^ ~ 3.85 



FIG. 5: Scaling function 6per{x) of the Casimir force for 
the d = 3 Heisenberg model in a slab geometry for peri- 
odic boundary conditions as function of the scaling variable 
X = tN^"^ , where A'^x ~ L±/a is the number of lattice lay- 
ers. The aspect ratio is chosen as N^^/N± = 6. Finite size 
scaling according to our expectation is confirmed within two 
standard deviations, where v = 0.71 has been chosen. The 
vertical scale has been adjusted according to the estimate 
^per(O) = 2Aper-,ri=3 = —0.78 (sBB main text). The error bars 
displayed here correspond to one standard deviation. 




FIG. 6: Variance Atx,x of the stress tensor for the Heisenberg 
model in d = 3 (see Eas. (l3.341 and 13.49ll 'l. normalized to 
A''^ at fixed aspect ratio N\^/N± = 6 as function of A'^x for 
different reduced temperatures t in the critical regime. The 
behavior is linear as indicated by the straight lines connecting 
the data points. Their slopes have been evaluated as 3.77 for 
t = -0.01, 3.92 for t = 0, and 4.03 for t = 0.01. The statistical 
error of the data (one standard deviation) is smaller than the 
symbol size. 



(see previous subsections and Eas. (|3.34|l and (|3.49|) ) for 
the spherical and the Gaussian model. The strict lin- 
earity also prevails for other temperatures in the scaling 
regime as shown in Fig^ In summary the Heisenberg 
model behaves just as the Ising and the XY model with 
respect to the variance of the stress tensor. 

Apart from different slopes there no specific differences 
in the behavior of the variance of the stress tensor for all 
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spin models investigated here in d = 3. However, the 
scahng function of the Casimir force does display specific 
differences as one may expect from the presence and the 
increasing dominance of Goldstone modes below Tc. 



V. SUMMARY AND CONCLUDING REMARKS 

In the current article an operator - the stress tensor op- 
erator - on a finite lattice systems has been constructed 
so, that its average value gives the universal behavior 
of the thermodynamic Casimir force near the critical 
point of a system with short-ranged interactions (see Eq. 
(|2.11(l 'l. The definition of the operator holds in systems in 
which the hyperscaling is valid (for 0{n) models that are 
systems with dimensionality 2 < d < 4) . Its explicit form 
for the two-dimensional Ising model is (see Eq. (|3.7() ') 



1 



(5.1) 



while that one for the c?-dimensional (2 < d < 4) spher- 
ical and the d-dimensional Gaussian models is (see Eqs. 
(pn^ and respectively) 



1 



'S'R'5'R+efc — (d— l)S'RS'R+ed 



.fc=l 



a V ^ 



(5.2) 



Here 7i is the Hamiltonian (normalized per unit particle) 
of the finite system, and Tij, is that one of the infinite sys- 
tem. For the spherical model one has to take into account 
that V — l/(d— 2), while z/ = 1/2 for the Gaussian model. 
On the example of the two-dimensional Ising model, the 
spherical model with 2 < d < 4 and the Gaussian model 
we have verified via exact calculations the correctness 
of the above presentation. They reproduce the correct 
values of the Casimir amplitudes at T — Tc and, for the 
spherical and the Gaussian models the expressions for the 
force derived via the excess free energy and via averag- 
ing the stress tensor operator are giving the same results. 
The amplitudes and the force near the critical point of 
the bulk system turns out, as expected, to be univer- 
sal and is in full accordance with the finite-size scaling 
theory. An evaluation of the variance of the so defined 
Casimir force has been also performed. If the summation 
is performed over all the particles within the system the 
corresponding result for the two-dimensional Ising model 
is (see Eq. (irTT|) ) 



Nj_Nn ~ 0.068Arj_Ar|, (5.3) 



while that one for the three-dimensional spherical and 
the Gaussian models is (see Eqs. (|3.34f) and H3.49|l ') 

At_L._L = A ^ t_L,_L(R-) - 0.107 A^j_iV|f . (5.4) 

R6A 

The average values of the above stress tensor operator 
are 

{T.t.Ahj))^-Y^ (5.5) 
for the two-dimensional Ising model (see Eq. (|3.4|l ). 

^ (E^^^^(i^))--^(^^^ii) (5-6) 



ReA 

0.306 , . 

=- {Nj_Nt 



(5.7) 



for the three-dimensional spherical model (see Eq. 
(ESU), and 



(U^^) ^ (^t^,^(R))=-l(|,(iV^Ar2) (5 



ReA 
0.382 

' Nl 



(5.9) 



for the three-dimensional Gaussian model (see Eq. 
|)). For the "noise-over-signal" ratio 



Pv 



(5.10) 



of the so-measured force from the above results one then 
derives 



PV ^ 0.159 ( ^ I N^, 



1/2 



for the Ising model, 



py ^ 1.069 ('^) 



for the spherical model, and 



PV c 0.856 ( ^ ) n'J' 



(5.11) 



(5.12) 



(5.13) 



for the Gaussian model. In the general case of a d- 
dimcnsional critical system the corresponding ratio at the 
bulk critical point is 



pv 



(d - 1)A V N\ 



(5.14) 



where D = D{T —t Tc) is a nonunivcrsal constant that 
describes the behavior of the variance of the tensor, i.e. 



(5.15) 
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FIG. 7; Scaling function 9per{x)/n of the Casimir force in 
d = 3 in a slab geometry for periodic boundary conditions as 
function of the scaling variable x = tN^'^ , where A'^j_ = L±/a 
is the number of lattice layers. Monte Carlo data are shown 
for the Ising model (+, n — 1), the XY model (x, n — 2) and 
the Heisenberg model (*, n — 3) with lattice size A^y = 180 
and N± = 30. The solid line shows 6per{x)/n in the spherical 
limit {n (x>). 



D[T) is a slowly varying nonuniversal function of T close 
to Tc, and A is the usual Casimir amplitude. 

Based on the proposed new operator, Monte Carlo cal- 
culations has been performed and the Casimir force scal- 
ing functions has been determined for the three dimen- 
sional Ising, XY and Heisenberg models. The scaling 
functions decay exponentially to zero above the critical 
temperature. The same happens for the Ising model also 
below Tc, while for the XY and Heisenberg models they 
tend to a constant because of the existence of the Gold- 
stone modes in this regime in these two models. Our 
results for 0{n) spin models, n = 1,2,3, oo, are summa- 
rized in Fig[7| The data for 9per{x) are normalized to n 
in order to obtain a direct comparison with the spherical 
limit, for which the exact result is shown. 

Our results confirm that one has to take into account 
the ratio between the thickness of the film and its lateral 
dimensions, when planning the settlement of an exper- 
iment, in order to achieve the desired noise-over-signal 
ratio. The numerical results that are presented can be 
considered as a type of such " measuring" of the force by 
Monte Carlo methods. They demonstrate clearly, that 
high accuracy in such type of measurement of the force 
is indeed possible to achieve. 
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APPENDIX A: THE CORRELATION LENGTH 
AMPLITUDE RELATION 

The coordinate transformation given by Ea. l|2.7|l re- 
moves the anisotropy from the spin model defined by 
Ea. (|2.1|) in the vicinity of the critical point. From the 
coordinate transformation and the principle of two scale 
- factor universality a relation between the correlation 
length amplitudes C||,o('^) £ind $_l,o(A) (see Ea. (|2.6|l ) can 
be established, which will be derived in the following un- 
der the assumption that hyperscaling is also valid, e.g., 
for 2 < c? < 4 and short-ranged interactions. 

According to the coordinate transformation given by 
Ea. (|^ we obtain ^[| = and = R{\) = (see 
£0. (12.81) 1. i.e., the parallel correlation length of the 
(untransformed) anisotropic system remains as the only 
correlation length of the (transformed) isotropic system. 
According to the principles of two scale - factor universal- 
ity and hyperscaling the singular part of the (bulk) free 
energy density ^^^^ [t) of the transformed spin system 
can then be written in the form 



-d 



(Al) 



where t = (3c{X)/(3 — 1 is the reduced temperature and 
A is a universal amplitude. Strictly speaking, one has 
to distinguish between two universal amplitudes A+ for 
T > Tc and A- for T < Tc. We disregard this distinction 
in Ea. (|Al|l in order to simplify the notation. Ea. (|Al|) is 
valid for T > Tc and T < Tc separately, provided, the 
correlation length remains finite for T < Tc- According 
to Ea. (|2.7l) the unit volume v of the system transforms 
as 



R{X)v 



(A2) 



and therefore we find 



f basing (-^7 ^) 



AR{\) [e||(A,t)]"' 



eii(A,t) 



= A[^^[\t)]-' [^11 (A, i)] 



-(rf-i) 



(A3) 



for the singular part of the bulk free energy density of 
the anisotropic, i.e., the untransformed system. 

According to Ea. l|2.6|l we have the alternative form 



fb,sing (A, i) 



-4[e±,o(A)]-^[eii.o(A)] 

A{X) \l3c{\)/p-l\'"' 



-id-i) 



\t\'^'' 
(A4) 



for Ea. ljA3l) . The nonuniversal amplitude ^(A) and /3c(A) 
must be independent of the labelling of the lattice axes, 
i.e., the direction which is chosen to be the 'perpendic- 
ular' one. From this symmetry argument and the par- 
ticular choice of the coupling constants J||(A) and J±{\) 
in Ea. H2.2|l we have already obtained /3'(A = 0) = in 
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Ea. (|2.5|) . Likewise, we obtain A'(A = 0) = from this 
symmetry argument. We therefore conclude that 



fb.sing (•^j t) 



= 



A=0 



and that due to 

A{X) ^ AiUA^T' [e||,o(Ar^'"'^ 
one also concludes from ^'(A = 0) that 



(A5) 



(A6) 



(A7) 



A=0 



From Eq. ljA7|l we finally obtain the important correlation 
length amplitude relation 



= (A8) 



A=0 



which is needed in the derivation of the stress tensor rep- 
resentation of the Casimir force for lattice spin models 
presented in Appendix IbI 



APPENDIX B: THE STRESS TENSOR 
REPRESENTATION OF THE CASIMIR FORCE 

The derivative of the excess free energy /ox with re- 
spect to A at the isotropic point A = is given by Eq. (|2.9() 
in the main text. The relation between Ea. H2.9|l and the 
Casimir force defined by Ea. (|l.l|) yields a lattice expres- 
sion of the stress tensor. This will be investigated here 
in the critical regime. Above the critical temperature all 
expressions will be exponentially small and can be ne- 
glected. Below the critical temperature Goldstone modes 
in 0{N > 2) systems also give rise to algebraically decay- 
ing finite-size effects, which will not be considered here. 

In order to find the relation between Eas. Hl.l|) and 
(|2.9I) we use the coordinate transformation given by 
Ea. (|2.7() and note that unlike the unit volume v (see 
Ea. (|A2|l ) the unit area remains invariant under Ea. H2.7|l . 
We recall that in the transformed (isotropic) system we 
have g(A) = g(A) = ^||,o('^) and we therefore find 

C(fex 

dX 



dX 



A=0 



A=0 



A=0 



dX 



df 



A-0 



_L,0 



dC, 



A=0 



dX 



c 



i?'(0) Lj 



A=0 

(Bl) 



_L.O 



A=0 



dX 



A=0 



where Eas. l|l.l) l and H2.8|l have been used. In order to 
evaluate the derivative df^^/d£,'j_ g in the critical regime, 
we use the critical finite-size scaling form 



(B2) 



and disregard the exponentially small contributions to 
Ea. (|Bl() from the regular part of the excess free en- 
ergy. Note, that for periodic boundary conditions the 
free energy of the finite system /'(t, L^) can be decom- 
posed, as usual, in a regular f^^g{t,L'^) and a singular 
/sing(^'^l) parts, where the regular part /^.^^{t, L'j_) can 
be taken to be equal (up to, eventually, exponentially 
small corrections) to that one of the infinite system, i.e. 
f regit: L'j_) — /rcg(i, oo) Q ■ That is why, for the periodic 
boundary conditions, the above equation ljB2|) is valid for 
the total excess free energy (and not only for its singular 
part). From Ea. l|B2p . we immediately obtain 



A=0 



ly £,0 dt ' 



(B3) 



where all terms on the r.h.s. of Eq.ljBSj) have already 
been evaluated at A = 0. To further evaluate Ea. l|B3|l we 
note that the excess internal energy Uox is given by 



5/c: 



/3c (0) du 



(B4) 



dl3 13^ dt 

In the vicinity oi (3 ~ /3c (0) Ea. ljB3|l can be rewritten as 



A=0 



1 



■(/3,(0)-/3) We 



(B5) 



In order to evaluate the derivative d£,'_^ q/ dX\\^Q in 
Ea. (|Bl(l we note that according to Ea. H2.7|l we have 
?lo(^) ~ C||,o(''^)- From the definition of R{X) given 
by Ea. H2.8|l we find by taking the derivative R'{X) with 
respect to A at A = 



1 










To 


dX 


A=0 


dX 


A=0- 



(B6) 



i?'(0) 



We eliminate j_,o/'^-^|a=o from Eg. lBSp using Ea. (|X^ 
of Appendix 1X1 and obtain 



^^11.0 



dX 



A=0 



I i?'(0). 



(B7) 



We finally insert Eas. (|B5l) and llBTl) into Ea. l|BT(l and, by 
rearranging the terms, we obtain for the Casimir force 



PFcas^■m^r = [i?'(0)] ' lim 



ill 



.k=l 



du L± 



(B8) 



where Ea. H2.9|l has also been used. Note that the first 
term in Ea. (jB8|) is generated by the anisotropy varia- 
tion whereas the second term originates from a change in 
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length scales enforced by the coordinate transformation 
given by Ea. (|^ . 

In order to express Ea. ljB8|l as the thermal average 
{t±±) of the normal component of the stress tensor we 
note that 



U — Ub 



(B9) 



where u is the volume energy density of the slab and ut 
the volume energy density in the bulk. Naturally, u and 
Ub are thermal averages of properly normalized Hamilto- 
nians H and Hb. More specifically, H is the Hamiltonian 
of the finite system normalized per unit volume, while 
Hb is the corresponding Hamiltonian for the bulk system 
(i.e. one imagines an arbitrary finite connected region of 
spins whose mutual probability distribution is obtained 
by taking the thermodynamic limit while integrating out 
all spins outside that fixed region. This is done for any 
finite region of the lattice). From Eas. l)B8|l and (|B9(I the 
operator form of the stress tensor given by Eq. (|2.11|) in 
the main text can then be read off. 



APPENDIX C: THE TWO-DIMENSIONAL ISING 
MODEL 

As it has been shown in the main text, see Eq. H3.7|l . 
in the critical region of the finite system the stress tensor 
is given by 

tx,x(i,j) — 2^/2^^^'''^^'''^^ ~ j-S'i+i j) 

+l{f3c-PKn^nb)- (ci) 

Let us now calculate the variance of the stress tensor 
Atx^x ih j): which we will interpret as a variance of a local 
measurement of the Casimir force made near the point 
(i,j). For the leading behavior of the variance near Tc 
one has 

= l-(5,,,+i5,+i,,)-(t,,,(*,j))'. (C2) 

Obviously, it holds that (Sij+iSi+ij) = (5'o,iS'i.o) = 
(50.05*1,1), because of the symmetry of the Ising model 
on a square lattice under periodic boundary conditions. 
The correlations (5o,o5i,i) are well known for the bulk 
system |0|: 
i) for T <Tc 



(50,05*1 1) — — E ( — 
TT \u 



(C3) 



ii) for r > Tc 



(5o,o^i,i) = — [E(u) + (u^ - 1)K(m)] . (C4) 



where, according to Eq. (|2.4|l . 



u = sinh(2/3J:r) sinh(2/3 J^). 



(C5) 



iii) for T = which is given by u = I, it follows that 
(5o,o5i,i) = 2/7r. 

In the above expressions K and E are the complete 
elliptic integrals of first, and of second kind, respectively. 
From them and Eq. ljC2|l one easily obtains expressions 
for the behavior of the variance ISix^xij-Tj) of the stress 
tensor below, above, and at T^. At T = Tc, for example, 
one has that 



Ai:,,^(z,j)~l-2/7r. 



(C6) 



Definitely, in addition from the above nonuniversal part 
the variance contains also an universal parts that are neg- 
ligible in comparison with the nonuniversal one. 

An estimation can be also derived for A ^ ■ ^. i3;,:c(i, j). 
With a variance of such a type one deals when, say, Monte 
Carlo simulations of the force are performed. According 
to Eqs. (|TT^ and (|TT3|l . at T ^ Tc 



A^L.x{iJ) 



AIL 



-/3n{\) 



(C7) 



^1 



-/3/(rc,A)], 



and, therefore, from Eq. H3.7|l it follows that 

I 



2 In (I 



71)"'"^' 



/3/(Tc,A)]. 



(C8) 

It is clear that the leading order behavior of the variance 
will stem from the bulk contribution to the free energy 
- the finite-size terms will produce only corrections to it. 
The bulk free energy of the anisotro pic two-dimensional 
Ising model is well known (see, e.g., |l5| ') 

-Pf ^ \n2+\f f If ln[cosh(2/3J.) 



X cosh{2j3Jy) — smh{2l3Jx) cos(6'i 
-sinh(2/3Jy)cos(6'2)] . 



(C9) 



Setting here = (1 + A) J and Jo, = (I — A) J, we im- 
mediately obtain —(3f{(3J,X), and from (|C8p one then 
derives (at T = T^) that 



(CIO) 



sin^(f )sin^(f ) 

-1 _ COs(fli)+COs(fl2) 

^ ^ 2 



1 2 

-2 + ^ 



ML ~ 0.068A/i. 



We again observe that the variance of the sum of tx^xihj) 
is proportional to the total number of summands in this 
sum. This is the result given in Eq. (|3.11|) in the main 
text. 
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APPENDIX D: THE SPHERICAL MODEL 

Using the identity 

ln(l + z)= / ^(1-e— )e-, 
Jo ^ 

the equation for the free energy 1)^11 7|l becomes 
where 



(Dl) 









+ sup < 


271 


w>0 L 



2N 

1 dx 

2 Jo ~ 



J(0) 



AT 



J(k)/J(0)] 



The supremum is attained at the value of w that is a 
solution of the (spherical field) equation 

^E /"e-^"'e-^[i-^('')/-^Wld:r = i^. (D4) 



keHA 



For nearest neighbor interactions the Fourier transform 
of the interaction reads 



J(k) ~ Jj cos{kjaj). 



(D5) 



Then, for the spherical field equation and the sum 
U{w,'N), we obtain 



TV. 



K, (D6) 



and 

U{w,N) = 
1 dx 



(DIO) 



d 

n 



qj—~oc 



In analogical way one can consider the behavior of the 
bulk system. Then, in the limit Nj ^ oo, j = I, d, 
one obtains the bulk equation for the spherical field 



(D2) K 
(D3) 



i=i 



(Dll) 



1 



Jo 



dni 



dUfi 



1 



w + S j=i ( 1 ~ Tij ) 

and the following contribution into the free energy 



1 dx 

2 Jo ~ 

1 1 



e--e--^Y[[e--''^Ioibjx)] 

3 = 1 



2 i^Trr Jo 
X In 



dni ■ 



dud 







w + ^ — cos TZj) 



(D12) 



In a similar way, starting from Eq. (|3.18|) . one can 
show that the bulk two-point correlation function in such 
an anisotropic system is 

1 ^ / h \ 

wPJ{0) Jo \ wj 



2^ 7o 



exp 



zn,Zo + p— cos n," 



.(D13) 



Supposing that Ij :$> 1, j — I,-- - ,d, from (|D13|I one 
obtains 



U{w,N) = 
1 f°° dx 



m G(r,<) 



1 



— cos k-i 



j=i ^ k, 



wPJ{0) Jo 



dp e-^n 



exp 



_w_j2 



1 VPbj/l 



, (D14) 



' wherefrom one concludes that the correlation length 
in direction j is 



where bj — Jj/J2'j=i Jj- Using the identity [Tsf 
2?™ 



(D15) 



N-l 



E ^'^p 



a; cos ■ 



TV 



= iv ^ V(a 



(D8) 



q— — oc 



Eqs. (jD6|) and (jD7|l can be written in the form 



-xw 



en 



with the critical point of the system given by w = 
(note that in the spherical model, because of the so-called 
equation of the spherical field, Eq. (jPlip . w depends on 
the coupling K, dimensionality d and on the anisotropy 
described by the constants hj, j — 1, ■ • ■ , c?). Therefore, 
one has 



dx = K, (D9) 



(Die) 
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Taking Jj, j = 1, • ■ ■ ,j in the form prescribed by Eq. 
(12. 2|) one obtains 



dX U±,o 



_ d J'jl) 
A=o " 2 J(l) ' 



(D17) 



and, thus, making use of Eq. (|2.11|) . we derive the explicit 
form of the stress tensor within the spherical model 



t±±(R) 



(3J 
1 



(D18) 



.fc=i 



a V ^ 



where 7i is the Hamiltonian (normalized per unit parti- 
cle) of the finite, and Ti;, of the infinite system. 



a. Evaluation of the finite-size excess free energy of the 
anisotropic system 



From Eqs. (|D10|I and (|D12(I one has 

U{w, N) = Ubiw) + AUiw, N) 

where 

d 



(D19) 



1 r°° lir 

^ ,_r. Jo 

(D20) 



q#0-" j=l 

Next, with the help of the expansion 

1 9 - 32i/2 



= exp(.-^/2x)^^^ 



^2ttx ^ 8a; 2!(8x)2 
AJ7(w, N) can be cast in the form 



) (D21) 



q#0 



,1 \/27rx6j 

(D22) 

wherefrom, in the limit of a film geometry 
iVi, N2, ■ ■ ■ Nd-i 00, with Nd — N±, one obtains 



Af7(u;,iVj_) = - 



(d-l)/2 



rf/4 



9=1 



(D23) 



Here we have taken Ji = J2 — ■ ■ ■ Ja-i = J\\ and Jd = 
J±, which corresponds to 61 = 62 = •■■^d-i = b\\ and 
bd — b±, whereas y = 2wN]^/b±. 

All what remains now is to deal with the behavior of 
the bulk term Ub{w) when K is close to Kc, i.e. when 
w << 1. This analysis is well known for the isotropic 
case, here we will, very briefly, extend it to cover the 



anisotropic case also. Starting from Eq. IID12|l . one ob- 
tains 



Ubiw) = UbiO) + 1 f dijWdHh), 
^ Jo 



where 







'loixbj) 



(D24) 



(D25) 



is a temperature independent constant and 

d 

Wd{Lu\h) = 



b) = / dxe-"'"r[e-"^^/o(a;6j) 
^0 j-^i 



(D26) 



is the generalized Watson-type integral (the standard one 
is with bj = 6 for all j = 1, • • • , d). Using the standard 
technique for evaluation of such type of integrals (see, 
e-g-, 0) one derives that, for 2 < c? < 4, 

WdHh) ~ Wd{0\h), (D27) 

wherefrom it follows that, again for 2 < d < 4, 



1 



Ubiw) ~ UbiO) + -wWd{0\h) 



1 



n~d/2) 



d/2 



2i2^)d/2UUVb] 

(D28) 

Taking into account that in terms of the "anisotropic" 
Watson integral the equation of the spherical field simply 
is 



K = Wdiw\h) 



(D29) 



and that, according to Eq. ljD15|l . the critical point is 
fixed by ui = 0, we conclude that the critical coupling of 
the anisotropic system is 

't' />oo d 

K, = 2(3,y2jj^WdiO\h)= / dxT\e-^''^Ioixbj). 

(D30) 

Then, close to if = Kc, for the free energy density of bulk 
system from Eqs. (ID2|) . HD28|I and H3.23|l one obtains 



PfbiK\h) 



2-K 



+ -WbiKc-K) + UbiO) 



1 



n-d/2) 



H27r)^/^UUVb: 



d/2 



(D31) 



where, for K < Kc, the parameter Wb is the solution of 
the equation 



K = Kc 



r(i - d/2) d/2-i 



(D32) 



whereas for K > Kc the supremum of the free energy is 
attained at it;;, = 0. Similarly, for the free energy density 
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of the finite system from Eqs. (|D2ll . (|D22|) . (|D28|) and 
H3.23|l we obtain 



Pf{K,N\h) = 



K 

In K 

27r 



+ -w{Kc -K) + Ub{0) 



1 



.,d/2 



(D33) 



where w is the solution of the finite-size equation for the 
spherical field 



K 



r(i -d/2) 



rf/2-l 



-J[ ^JQjtxFj 



poo d 
q#0-^° 3 = 1 



(D34) 



Recalling that, for 2 < d < 4, the spherical model has a 
critical exponent v = l/('^ — 2) one can, in the limit of a 
film geometry A^i, N2^ • • • , Nd~i — > oo, from Eqs. l|D23p . 
(|D3ip and l|D33|) , obtain an expression for the excess free 
energy f3{f - ft) (per spin) 



(D35) 



ir(-d/2) fb^\'''^'^^^^ 



2 (47r)'^/2 
(27r)'i/2 V b 



(^) (^^^^-^^^) 

y 1. j 



in a scaling form. In the above equation 



2' 



(D36) 



is the temperature scaling variable, yoa = 2wbN'j_/b± is 
the solution of the bulk spherical field equation 



1^ _ r(l-d/2) fbj_ 



(47r)'^/2 



(d-l)/2 



yd/2-1 



(D37) 



while y — 2wN'^/b± is the solution of the finite-size 
spherical field equation 



{47rY/2 \ 5|| J 



(d-l)/2 



y 



d/2-l 



(D38) 



(27r)'^/2 \ b 



(d-l)/2 



d/4-1/2-^ ^d/2-l(gv^) 



9=1 



ci/2-1 



For the Casimir force (see also Eq. (11.11) 1 
9 



[N±Kf - A)] 



(D39) 



from Eqs. (jX^ . l|X77jl one obtains 



c asimir 



^I''<i ^2^1(2/ -yoo) 



(D40) 



-id-l) 



(27r)<^/2 



(d-l)/2 



1 r(-rf/2) ^^^,/2 

2 (47r)''/ 



./2„ 

^W/2 2/00 J 



d/4 \ " 



E 

9=1 



jrf/2 



We are ready now to determine the Casmir amplitudes 
A in the spherical model. Having in mind Eq. (|1.4|l . at 
K = Kc (then j/oo = 0), for the isotropic system (then 
6_L = b\\ = 1/d) one obtains from Eq. I|3.28(l 



A 



where 



ir(-d/2) 
2 (47r)rf/2 ' 



d/2 



(27r)''/2 



d/4 
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■7=1 



Kd/2{q^/y^) 



rjd/2 



(D41) 
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(47r)''/2 d(27r)''/2^= 
Using now that (see, e.g., [1^ 1 



V2^. 

9=1 



Kd/2-i{qJv'c) 



rjd/2-l 



2u 



K,+i{z)^K,^i{z) + —K,{z), 
from Eqs. l|D4n) and ||D42|) we derive 



A = 



d(27r)^/2 



4+1/2 



Kd/2+iiq^/y^) 



q=l 



q' 



d/2-1 



(D42) 
(D43) 

(D44) 



with /3cf'casimir(ifc,i±) = (d " 1) A Lj"^. Thc cxact 
value of ?/c and A in an explicit form is only known for 
d = 3. Then 



(D45) 



y = 41n"[(H- V5)/2 



(this value is well-known and seems that has been derived 
for the first time in 20] ), and pll | 



A 



2C(3) 



(D46) 



This is the only exactly know Casimir amplitude for a 
three dimensional system. In |2^ it has been shown (see 
there Eq. (27)) that this value can also be written in the 
form 



A. -A 
2tt 



D 



(D47) 

where Lip(z) = J2T=i z'^/k^ is the polylogarithm func- 
tion of order p. Taking into account that K{5/2,x) — 
^7r/2a;(l -|- 2/x + 3/x'^) exp{-x) and Eq. one 
can easily check that the right hand side of Eq. 
can indeed be written in the form given in (jD47|l 
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b. Evaluation of the average value of the stress tensor 

Taking into account H2.16|l . for the difference of the 
finite-size and bulk internal energy densities u = (7i) 
and lift = (Tif,), one can easily derive from ljD2|) 



U — Ub 



-J{W - Wb) 



1 

Ad 



J{y-yoo)N-\, (D48) 



where y = 2dwNl, and then from (|3.18|) . or (jD2p . to 
obtain for the stress tensor (13.221) that 



1 1 



d N 



^-v cos(/ciai) — cos(fcrfac;) 
keBA + - Ej=i cos(fcjaj) 



rf- 2 



(D49) 



4d 

where, we recall, xi = d~^{Kc — K)N'[~'^ (see Eq. 



Using the identity (see Eq. (|D8 
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\ N J 

(D50) 

where /^.(a;) = '3l^-^v{x), in the limit of a film geometry, 
i.e. when Ni, N2, - ■ ■ , oo, the above expression 

can be rewritten in the form 



{t±±m = 



2(d- 1) 
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-2x 
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[e-^/o(a 
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id 

It is worth to mention that till now no approximation in 
the calculation of the of the average of the stress tensor 
operator has been made. In order to obtain the scaling 
form of the above expression such a step will be per- 
formed only now. Indeed, with the help of the expansion 
ID2T1 



) (D52) 



I^{x) = ' ^(1 + — - 



/2Trx ' 8x 2l{8x)^ 
one can set the above expression in scaling form 



rf(27r)'i/2 
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(D53) 



Now it only remains to show that the right-hand side of 
the above equation is indeed equal to the right-hand side 
of Eq. H3.28|l (for b± = 5||), which gives the Casimir force 
calculated in a direct manner as a derivative of the finite- 
size free energy with respect of the size of the system. In 
order to demonstrate that, let us first, with the help of 



identity ljD43|) , rewrite the above expression for the stress 
tensor in the form 



2{d~l) 

{27Ty/^ 
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Next, from the bulk (|3.2()|) and finite-size equations ()3.27fl 
of the spherical field one directly derives 
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-^xiyoc 
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wherefrom 



(27r)'i/2 "" " E 'qd/2-l ' - 4^1(2^ 



Inserting now ljD57|l in (|3.32|l we conclude, that, indeed. 



a 

for the spherical model. 



(D58) 



c. Evaluation of the variance of the stress tensor 

If AC is the variance of the random variable C, i.e. 
AC =< (C- < C >f >=< > - < C >^ then at 

T = n 



&[^"E' 



-l3n(X) 



(D59) 



(for the definition of n{X)) and /(T, A) see Eqs. (ITT^ - 
(|2.15(l . Comparing now Eqs. H2.15|l and (|3.22(l . we con- 
clude that, at T = Tc, in the case of a spherical model 



A ^ f^,^(R) = ±N^N^-'^[-(3fin,X)]. (D60) 
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The finite-size free energy density of the anisotropic sys- 
tem is given in Eq. (|D33|I . where the anisotropy is char- 
acterized by the constants 



wherefrom one immediately derives 



foi = &2 = • • • = bd-i = &II 
and 

bd = b± 



Jii 



1 + A 



(d- 1)J|| + Jj_ 



(d- 1)J|| + Jj_ 
Let us now first note that 



d-1 
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Performing now the calculations, from l|D25|l and (^1 
we obtain 



q2 ^ roo 

K = /3J(0) = 2/3 [{d - 1) Jji + J ±] = 2/3Jd (D63) ^t^t'(0|b) = -d{d - 1) dxxe-'''-' I^^^x) (D68) 



does not depend on A and that 
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It is clear from Eqs. HD33|) and (|3.24|l ') that the con- 
tributions to the variance of the stress tensor stemming 
from the "finite size" and the singular "bulk" parts will 
be of the order of {N±N^~'^)/Nf, wile that one from 

the "bulk" regular part will be of the order of Nj_N^^^. 
Because of that the leading contributions will be nonuni- 
versal. In addition to them one will have also universal 
corrections, but we will neglect them in the current treat- 
ment and will deal only with the leading-order behavior 
of the variance of the Casimir force. Then, from Eq. 
(|D33|I . one has 



which leads to the following result 



ReA 



dxxlg-^ix) 



-/o(x)(/o(x)+/2(a;))^e-'^^ (D69) 
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-f/b(0|b], (D65) 



for the variance of the Casimir force in the spherical 
model at r = Tc. This will be also the leading result 
everywhere in the critical region. A numerical evalua- 
tion of Eq. HD69|I gives 



where £/f,(0|b ) = £ /fc(0) is defined in Eq. |D25)l . Having 
in mind Eq. ljD64|) , it is easy to show that 



AU,± = A ^ U,_l(R) ^ 0.107 Nx^Nl 



(D70) 
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